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Abstract 

The study of the ultrarelativistic plasmas in perturbation theory is plagued with 
infrared divergences which are not eliminated by the screening corrections. They affect, 
in particular, the computation of the lifetime of the elementary excitations, thus casting 
doubt on the validity of the quasiparticle picture. We show that, for Abelian plasmas 
at least, the infrared problem of the damping rate can be solved by a non-perturbative 
treatment based on the Bloch-Nordsieck approximation. The resulting expression of 
the fcrmion propagator is free of divergences, and exhibits a non- exponential damping 
at large times: Sij(i) ~ exp{— aTt lnw p £}, where oj p — gT/3 is the plasma frequency 
and a = g 2 /Air. 
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The study of the uitrareiativistic plasmas in perturbation theory is plagued with infrared divergences 
which are not eliminated by the screening corrections. They affect, in particular, the computation 
of the lifetime of the elementary excitations, thus casting doubt on the validity of the quasiparticle 
picture. We show that, for Abelian plasmas at least, the infrared problem of the damping rate 
can be solved by a non-perturbative treatment based on the Bloch-Nordsieck approximation. The 
resulting expression of the fermion propagator is free of divergences, and exhibits a non-exponential 
damping at large times: Sn(t) ~ exp{— aTt \nu} p t}, where ui p = gT/3 is the plasma frequency and 
a — g /4-7T. 



1 Introduction 

The study of the elementary excitations of uitrareiativistic plasmas, such as the quark- 
gluon plasma, has received much attention in the recent past u'EI. The physical picture 
which emerges, for both Abelian and non-Abelian gauge theories, is that of a system with 
two types of degrees of freedom: i) the plasma quasiparticles, whose energy is of the order 
of the temperature T; ii) the collective excitations, whose typical energy is gT, where g is 
the gauge coupling, assumed to be small: g <C 1 (in QED, g = e is the electric charge). 

For this picture to make sense, however, it is important that the lifetime of the excita- 
tions be large compared to the typical period of the modes. Information about the lifetime 
is obtained from the retarded propagator. A usual expectation is that Sn(t, p) decays expo- 
nentially in time, Sn(t, p) ~ e - lE (p) t e -i{p) t ^ where E(p) ~ T or gT is the average energy of 
the excitation, and jip) is the damping rate. The exponential decay may then be associated 
to a pole of the Fourier transform Sr(uj, p), located at lo = E — £7 : 

/oo 1 
dte- wt S R (t,p) (1) 
-00 w - {E{p) - i7(p)) 

The quasiparticles are well defined if the damping rate 7 is small compared to the energy 
E. If this is the case, the respective damping rates can be computed from the imaginary 
part of the on-shell self-energy, £(u> = E(p),p). 

Previous calculations a suggest that 7 ~ g 2 T for both the single-particle and the collec- 
tive excitations. In the weak coupling regime g <C 1, this is indeed small compared to the 
corresponding energies (of order T and gT, respectively), suggesting that the quasiparticles 
are well defined, and the collective modes are weakly damped. However, the computation 
of 7 in perturbation theory is plagued with infrared (IR) divergences, which casts doubt on 
the validity of these statements 3-7 . 

The occurence of IR divergences in the perturbation theory for hot gauge theories is 
generic, and reflects the long-range character of the gauge interactions. However, it is well 

''Laboratoire de la Direction des Sciences de la Matiere du Commissariat a l'Energie Atomique 



known that most of these divergences are actually removed by screening effects which are 
generated by the collective motion of the thermal particles. In ultrarelativistic plasmas, the 
screening effects manifest themselves over typical space-time scales ~ 1/gT. Their inclusion 
in the perturbative expansion — which is achieved in a gauge-invariant way by resumming 
the so called "hard thermal loops" (HTL) of Braaten and PisarskiH — greatly improve the 
infrared behavior, and yields IR- finite results for the transport cross-sections oa. and also 
for the damping rates of the excitations with zero momentum tH12l. 

However, the HTL resummation is not sufficient to render finite the damping rates of 
the excitations with non vanishing momentatl The remaining infrared divergences are due 
to collisions involving the exchange of longwavelength, quasistatic, magnetic photons (or 
gluons), which are not screened in the hard thermal loop approximation. Such divergences 
affect the computation of the damping rates of charged excitations (fermions and gluons), in 
both Abelian and non-Abelian gauge theories. Furthermore, the problem appears for both 
soft (p ~ gT) and hard (p ~ T) quasiparticles. In QCD this problem is generally avoided 
by the ad-hoc introduction of an IR cut-off ("magnetic screening mass") ~ g 2 T, which is 
expected to appear dynamically from gluon self- interactions In QED, on the other hand, 
it is known that no magnetic screening can occur 0, so that the solution of the problem 
must lie somewhere else. 

We have shown recently^ that, for Abelian plasmas, the divergences can be eliminated 
through a non perturbative treatment, which involves a soft photon resummation a la 
Bloch-Nordsiecklij. We have thus obtained the large-time decay of the fermion propagator, 
which is not of the exponential type alluded to before, but of the more complicated form 
Sil(t) ~ e ~ lEt exp{— aTt Inupt}, where a = g 2 /^ and ui p ~ gT is the plasma frequency. 
Accordingly, the Fourier transform Sr(u) is analytic in the vicinity of the mass-shell. In 
what follows, I will briefly discuss these results, their derivation and their consequences. 



2 Perturbation theory for 7 

Let me first recall how the infrared problem occurs in the perturbative calculation of the 
damping rate 7. For simplicity, I consider an Abelian plasma, as described by QED, and 
compute the damping rate of a hard electron, with momentum p ~ T and energy E(p) = p. 

To leading order in g, and after the resummation of the screening corrections, 7 is ob- 
tained from the imaginary part of the effective one-loop self-energy in Fig. |. The blob 
on the photon line in this figure denotes the effective photon propagator in the HTL ap- 
proximation, commonly denoted as *D fJjU (q). In the Coulomb gauge, the only non-trivial 
components of *D fMU (q) are the electric (or longitudinal) one *Z?oo( ( z) = *Di{q), and the 
magnetic (or transverse) one *Dij(q) = (5ij — qiqj)*D t (q), with 

*Di(q ,q) = — ^ r, *D t (q ,q) = — 2 — \— -, (2) 

q 2 -Ui{q ,q) q£ - q 2 -n t {q ,q) 

where and Tit are the respective pieces of the photon polarisation tensor S'i. Physically, 
the on-shell discontinuity of the diagram in Fig. |] accounts for the scattering of the incoming 
electron off the thermal fermions, with the exchange of a soft, dressed, virtual photon. 



p 
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Figure 1: The resummed one- loop self-energy 



The corresponding interaction rate is simply computed as 7 = op, where p ~ T 3 is the 
density of the scatterers, and a = J d 2 q (da/dq 2 ), with q denoting the momentum of the 
exchanged (virtual) photon. For a bare (i.e., unscreened) photon, the Rutherford formula 
yields da/dq 2 ~ g 4 /q 4 , so that 7 ~ g 4 T 3 J(dq/q 3 ) is quadratically infrared divergent. 
Actually, the screening effects included in IIj,t soften this IR behaviour. We have, in the 
kinematical regime of interest, (uo p = eT/3) 

M% « Q) - H = m 2 D , U t (q « q) ~ -i oo 2 p ^ . (3) 

We see that screening occurs in different ways in the electric and the magnetic sectors. In 
the electric sector, the familiar static Debye screening provides an IR cut-off mn ~ gT. 
Accordingly, the electric contribution to 7 is finite, and of the order 7/ ~ g 4 T 3 /m 2 D ~ g 2 T. 
In the magnetic sector, screening occurs only for nonzero frequency q^. This comes from 
the imaginary part of the polarisation tensor, and corresponds to the absorption of the 
space-like photons < q 2 ) by the hard thermal fermions (Landau damping ta). This 
"dynamical screening" is not sufficient to completely remove the IR divergence of 7t , which 
is just reduced to a logarithmic one: 



[■00 rq 

lt ~ g 4 T 3 / dq / dq \*D t (q ,q)\ 2 
JO J-g 



00 rq 



>v dq 



With an IR cut-off p, "ft ~ g 2 T^ n {'^p/ p)- The remaining logarithmic divergence is due to 
collisions involving the exchange of very soft, quasistatic (qo — > 0), magnetic photons, which 
are not screened by plasma effects. To see that, note that the IR contribution to 7t comes 
from momenta q <C gT, where \*Df(qo, q)\ 2 is almost a delta function of q^: 

\*D t (q ,q)\ 2 ~ 4 1 2 — 7-^2 > q ^o iT^S(qo). (5) 

This is so because, as qo — > 0, the imaginary part of the polarisation tensor vanishes linearly 
(see the second equation (||)), a property which can be related to the behaviour of the phase 
space for the Landau damping processes. Since energy conservation requires qo = qcos8, 
where 6 is the angle between the momentum of the virtual photon (q) and that of the 
incoming fermion (p), the magnetic photons which are responsible for the singularity are 
emitted, or absorbed, at nearly 90 degrees. 




Figure 2: A generic n-loop diagram (here, n — 6) for the self-energy in quenched QED. 
3 A non perturbative calculation 

The IR divergence of the leading order calculation invites to a more thorough investigation 
of the higher orders contributions to 7. Such an analysis Bill reveals strong, power-like, 
infrared divergences, which signal the breakdown of the perturbation theory. To a given 
order in the loop expansion, the most singular contributions to 7 arise from the quenched 
(no internal fermion loops) self-energy diagrams in Fig. 0, where all the internal photon 
lines are of the magnetic type. As in the one loop calculation, the leading divergences arise, 
in all orders, from the kinematical regime where the internal photons are soft (q — > 0) and 
quasistatic (qo — ► 0). Physically, these divergences come from multiple magnetic collisions. 

This peculiar kinematical regime can be conveniently exploited in the imaginary time 
formalism @, where the internal photon lines carry only discrete (and purely imaginary) en- 
ergies, of the form qo = ioo n = i2irnT, with integer n (the so-called Matsubara frequencies). 
Note that the non-static modes, with n / 0, are well separated from the static one qo = 
by a gap of order T. In this formalism, all the leading IR divergences of the damping rate 
- which, I recall, arise from the kinematical limit qo — > — are concentrated in diagrams 
in which the photon lines are static, i.e., they carry zero Matsubara frequency ill!. Note 
that, for these diagrams, all the loop integrations are three-dimensional (they run over the 
three-momenta of the internal photons), so that the associated IR divergences are those of 
a three-dimensional gauge theory. This clearly emphasizes the non perturbative character 
of the leading IR structure. 

In what follows, we restrict ourselves to these diagrams, and compute their contribution 
to the fermion propagator near the mass-shell, in anon perturbative way. This can be 
"exactly" done in the Bloch-Nordsieck approximationta, which is the relevant approximation 
for the infrared structure of interest. Namely, since the incoming fermion is interacting only 
with very soft (q — > 0) static (qo = 0) magnetic photons, its trajectory is not significantly 
deviated by the successive collisions, and its spin state does not change. Thus, we can ignore 
the spin degrees of freedom, which play no dynamical role, and we can assume the fermion 
to move along a straightline trajectory with constant velocity v (for the ultrarelativistic 
hard fermion, |v| = 1; more generally, for the soft excitations, v(p) = dE(p)/dp = v(p)p is 
the corresponding group velocity, with \v(p)\ < 1). Under these assumptions, the fermion 



propagator can be easily computed as@ 

S R (t,p) = i9(t)e- iE ^A(t), (6) 

where 

Am / 2^ r p d 3 g 1 l-cost(v(p) -q) | 

contains all the non-trivial time dependence. The integral in eq. (j7|) is formally identical 
to that one would get in the Bloch-Nordsieck model in 3 dimensions. Note, however, the 
upper cut-off u> p ~ gT, which occurs for the same reasons as in eq. (Q). Namely, it reflects 
the dynamical cut-off at momenta ~ gT, as provided by the Landau damping. 

The integral over q has no infrared divergence, but one can verify that the expansion of 
A(i) in powers of g 2 generates the most singular pieces of the usual perturbative expansion 
for the self-energy El. Because our approximations preserve only the leading infrared be- 
havior of the perturbation theory, eq. ([?]) describes only the leading large-time (t S> 1/gT) 
behavior of A(t). This is gauge independently and of the form (we set here a = <? 2 /47r and 
v(p) = 1 to simplify writing) 

A(uj p t > 1) ~ exp(-aTt]nu p tY (8) 

Thus, contrary to what perturbation theory predicts, A(i) is decreasing faster than any 
exponential. It follows that the Fourier transform 

/oo /*oo 
dte-^S R (t,p) = i \ dte lt ^~ E ^ + ^A(t), (9) 
-oo JO 

exists for any complex (and finite) u. Thus, the retarded propagator Sr(uj) is an entire 
function, with sole singularity at Imu — > — oo. The associated spectral density p(ui,p) 
(proportional to the imaginary part of Sr(u,p)) retains the shape of a resonance strongly 
peaked around the perturbative mass-shell u = E(p), with a typical width of order ~ 
g 2 T\n(l/g)B. 



4 Conclusions 



The previous analysis solves the IR problem of the damping rate, thus confirming the 
quasiparticle picture for the hot Abelian plasmas. For high temperature QCD, on the other 
hand, the resolution of the corresponding problem also requires the understanding of the 
non-perturbative sector of the magnetostatic gluons. 

As a crude model of QCD, let us assume, in agreement with heuristic argumentsB, and 
also with lattice computations @, that a screening mass [i ~ g 2 T is dynamically generated 
in the magnetic sector. After replacing 1/q 2 — > l/(q 2 + fJ- 2 ) for the photon propagator 
in eq. (|7|), this equation provides, at very large times t ^ l/g 2 T, an exponential decay: 
A(i) ~ exp(— jt), with 7 = aT\n(cj p / fi) = aT\n(\/g). However, in the physically more 
interesting regime of intermediate times 1/gT < t < l/g 2 T, the behavior is non-exponential 



and governed by the plasma frequency, according to our result (|8|): A(i) ~ exp(— aTt lnw^t). 
Thus, at least within this limited model, which is QED with a "magnetic mass", the time 
behavior in the physical regime remains controlled by the Bloch-Nordsieck mechanism. But, 
of course, this result gives no serious indication about the real situation in QCD, since it is 
unknown whether, in the present problem, the effects of the gluon self-interactions can be 
simply summarized in terms of a magnetic mass. 
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